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statistical  techniques  are  used,  time/frequency  domain  analysis  of  the  IF 
signal  is  not  required  and  FM/FM  systems  or  systems  which  use  random  modulations 
are  conveniently  analyzed.  The  range  response  is  derived  in  terms  of  a sum  of 
weighted  quasi-autocorrelation  functions  which  have  been  translated  along  the 
delay  axis.  The  displacements  and  weightings  depend  on  the  harmonic  series 
which  describes  the  IF  reference  as  a function  of  the  frequency  modulation. 

The  analysis  technique  is  illustrated  with  examples  of  systems  designed  for  - 
low-height  ranging  applications.  Calculated  and  measured  results  are  shown  for 
simple  hybrid  modulation/demodulation  systems  synthesized  using  the  techniques 
described  in  the  report.  These  systems  feature  power  spectra  which  yield  well 
defined  range  responses  and  may  employ  FM/FM  techniques  to  reduce  the  frequency 
ambiguities  inherent  in  harmonic  systems. 
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STATISTICAL  ANALYSIS  OF  A CLASS 
OF  IF  CORRELATOR  FM  RANGING  SYSTEMS 


1.  INTRODUCTION 


This  report  describes  a statistical  technique  for  analyzing  a class  of  FM 
ranging  systems  called  "IF  Correlators."  IF  correlators  are  defined  as  systems 
which  correlate  the  IF  signal  which  results  from  a target  at  delay  time  t with 
a locally  generated  reference  which  is  an  estimate  to  the  IF  signal  expected  at 
t = tr,  where  rr  is  a reference  delay.  The  range  response  of  IF  correlators 
maximizes  near  t = tr  and  will  have  range  resolution  determined  by  the  transmitted 
RF  bandwidth  B.  The  IF  correlation  technique  has  been  considered  for  a variety  of 
ranging  applications.  References  [l-7]  analyze  various  aspects  of  these  systems. 


•i 

I, 


The  approach  generally  used  for  analyzing  IF  correlators  is  to  describe  the 
IF  signal  s(t)  and  reference  r(t)  as  deterministic  time  functions.  The  correlator  i 

low-pass  output  is  then  the  DC  term  resulting  from  expressing  the  product  s(t)*r(t)  l 

in  a Fourier  (time)  series.  Target  delay  t is  treated  as  a parameter  since  it  is 
slowly  varying;  the  variation  of  the  correlator  output  due  to  varying  t produces  f 

the  range  response.  For  simple  periodic  FM  modulations  the  time  analysis  approach  ) 

leads  readily  to  useful  results;  however,  this  approach  is  very  difficult  for  IF  | 

correlators  which  use  random  or  periodic  FM/FM  type  modulations  [8].  The  time  | 

analysis  also  obscures  important  functional  relationships  which  are  time-indepen- 
dent. 


The  statistical  analysis  technique  described  in  this  report  is  based  on  the  ,, 

fact  that  both  the  IF  signal  and  reference  can  be  expressed  as  sinusoidal  functions  J 

(or  a sum  of  sinusoidal  functions)  of  a common  angle  0 which  is  directly  related  | 

to  the  instantaneous  transmitted  frequency.  The  angle  0 is  the  phase  difference  ! 

between  the  transmitted  and  received  signals  evaluated  at  a particular  time  delay  t l 

and,  for  small  delays  and  slowly  varying  modulations,  is  simply  the  product  of 
the  instantaneous  radian  frequency  and  delay  time  t.  The  DC  value  resulting  from  < 

the  product  of  the  signal  S (G ) and  reference  R(0)  is  found  by  statistically  aver- 
aging the  product  using  the  probability  density  of  0,  p(0).  0 is  directly  related 

to  the  instantaneous  transmitted  frequency;  hence,  for  high-index  FM  modulation, 
p(0)  is  directly  proportional  to  the  envelope  of  the  transmitted  power  spectrum.  i 

Since  R(0)  and  S(0)  are  sinusoidal  (or  a sum  of  sinusoids),  their  product  is  sinus- 
oidal in  0.  Statistically  averaging  the  sinusoids  resulting  from  the  product 
S(0)*R(0)  using  p(0)  (which  is  related  to  the  power  spectrum)  is  equivalent  to 
Fourier  transforming  p(0).  As  a consequence,  the  resulting  range  response  is  ex- 
pressed in  terms  of  the  autocorrelation  of  the  transmitted  RF  signal. 

The  statistical  technique  depends  on  a time-independent  description  of  the 
signal  and  reference  as  a function  of  the  angle  0 and  is  therefore  applicable  to 
random,  FM/FM,  and  periodic  modulations  with  equal  facility.  Moreover,  it  provides 
a unifying  concept  for  analyzing  and  comparing  the  variety  of  systems  described  in 
the  cited  r: terences  by  demonstrating  that  all  IF  correlator  systems  which  have 
the  same  modulation/reference  statistics  will  have  the  same  range  responses.  This 
fact  is  by  no  means  apparent  using  time  domain  techniques  since  equivalent  systems 
may  appear  quite  different. 


The  statistical  approach  has  proven  useful  for  the  design  and  analysis  of  a 
variety  of  IF  correlator  FM  ranging  systems.  In  most  cases,  nearly  optimum  sys- 
tems can  be  designed  using  simple  graphical  techniques  and  their  range  responses 
expressed  in  closed  form.  This  technique  shows  that  all  IF  correlators  (harmonic 
N systems,  etc.)  are  directly  related  to  IF  replica  correlators  which  use  an  RF 
delay  line/mixer  to  produce  an  IF  replica.  (A  replica  is  an  exact  copy  of  the 
IF  signal  expected  from  an  ideal  target  at  t * ir.)  It  also  demonstrates  that 
any  IF  reference  can  be  expressed  in  terms  of  replicas  and  the  range  response 
expressed  in  terms  of  weighted  displaced  autocorrelation  functions.  The  weighting 
and  displacements  are  related  to  the  Fourier  series  which  defines  the  reference  R 
as  a function  of  8;  the  autocorrelation  function  is  found  by  transforming  p(Q). 

The  statistical  averaging  technique  used  is  quite  general  and  is  dependent 
only  on  an  assumption  of  ergodicity  (the  equivalence  of  time  and  statistical  aver- 
aging) which  is  satisfied  (approximately)  for  most  physical  systems.  To  simplify 
the  analyses  to  follow,  certain  assumptions  which  are  reasonable  for  the  class  of 
FM  ranging  systems  being  investigated  will  be  made.  First,  high-index  FM  modu- 
lation of  the  oscillator,  with  no  incidental  AM,  is  assumed.  For  this  case,  the 
envelope  of  the  transmitted  power  spectrum  will  be  directly  proportional  to  the 
probability  density  of  the  frequency  modulation.  Second,  a linear  voltage-to- 
frequency  FM  modulation  characteristic  is  assumed.  It  then  follows  that  the  pro- 
bability density  of  the  frequency  modulation,  and  hence  the  envelope  of  the  trans- 
mitted power  spectrum,  will  be  directly  proportional  to  the  probability  density 
of  the  modulating  voltage.  Third,  the  probability  density  of  the  modulation,  p(0), 
will  be  assumed  to  be  an  even  function,  implying  an  even  transmitted  power  spectrum 
about  the  center  frequency  u)0.  This  assumption  simplifies  the  mathematical  expres- 
sion for  the  autocorrelation  function.  Fourth,  only  ranging  systems  for  which  the 
product  of  the  IF  center  frequency  and  the  target  time  delay  tr  is  small  will  be 
considered.  With  this  assumption,  certain  phase  errors  can  be  ignored  and  the 
range  resj  'mse  analyzed  using  first-order  statistics. 

Using  the  above  assumptions,  the  signal  can  always  be  expressed  as  a single- 
valued sinusoidal  function  of  8.  For  IF  correlators,  the  reference  is  necessarily 
synchronized  to  the  modulation;  therefore,  the  functional  relationship  R(8)  is  ex- 
plicitly determined  by  the  hardware  implementation.  When  R(0)  is  not  sinusoidal 
in  0,  it  can  always  be  expressed  in  a Fourier  series  as  a weighted  sum  of  sinusoids. 
For  most  IF  correlators  R(0)  will  be  a singled-valued  function;  however,  for  some 
systems  it  will  be  multi-valued  as  will  be  illustrated  in  the  examples.  When  the 
reference  is  multi-valued,  the  modulation  is  partitioned  into  discrete  monotonic 
sections  for  which  the  reference  is  single-valued.  The  system  range  response  may 
then  be  computed  by  summing  the  range  responses  due  to  each  monotonic  section; 
alternatively,  a composite  R(0)  can  be  determined  by  averaging  the  multiple  values 
according  to  their  probabilities.  The  system  range  response  is  then  expressed  in 
terms  of  the  autocorrelation  function  obtained  from  the  Fourier  transform  of  the 
total  transmitted  power  spectrum. 

Section  2 describes  IF  correlator  FM  ranging  systems  and  presents  some  back- 
ground material.  Section  3 describes  the  IF  signal  and  reference  in  statistical 
terms  and  presents  the  necessary  justification  for  the  transition  from  time  to 
statistical  averaging.  Section  4 derives  the  range  response  using  the  statistical 
descriptions  of  Section  3.  Section  5 describes  a procedure  for  computing  the  range 
response  using  the  statistical  technique.  This  procedure  is  used  to  compute  the 
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range  response  for  numerous  examples  which  are  chosen  to  illustrate  the  major 
concepts  presented  in  the  report.  Section  6 presents  several  simple  hybrid 
modulation  systems  which  were  used  to  obtain  test  results  and  compares  computed 
and  measured  range  responses.  Section  7 presents  the  conclusions. 

2.  IF  CORRELATOR  FM  RANGING  SYSTEM  DESCRIPTION. 

Systems  are  described  by  the  block  diagram  of  Figure  1,  which  shows  a single 
antenna  configuration  and  envelope  detection  of  the  composite  RF  voltage  at  the 
antenna.  Normalized  signals,  references,  autocorrelations,  etc.,  will  be  used 
to  simplify  the  analysis.  The  output  of  practical  systems  will,  of  course,  vary 
with  signal  amplitudes,  mixer  gain  coefficients,  etc. 


Figure  1:  IF  Correlator  FM  Ranging  System 


The  IF  signal  s(t)  is  multiplied  by  a reference  r(t)  and  the  product  fil- 
tered (time-averaged)  to  obtain  the  range  response  RL(t).  RL(t)  can  be  written 
as  the  pioduct  of  an  amplitude  factor  usually  called  the  range  law  and  a sinusoidal 
doppler  carrier.  Variations  in  target  delay  t will  cause  variations  in  RL  at  dop- 
pler  frequencies  which  are  much  lower  than  the  modulation/demodulat inn-system  and 
RF-carrier  frequencies;  therefore,  r can  be  considered  a slowly  varying  parameter 
in  the  analysis  to  follow. 


Figure  1 illustrates  that  the  reference  may  be  related  to  the  modulation  ir 
two  ways.  First,  the  reference  may  be  derived  from  the  RF  signal  using  a technique 
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which  is  equivalent  to  the  RF  delay  line/mixer  described  by  Peperone  [l] . For 
this  case,  the  reference  will  have  information  about  oscillator  center  frequency 
u)q  as  well  as  the  varying  component  <%(t).  A reference  of  this  type  (replica)  is 
a proper  IF  reference  for  all  modulations  and  time  delays  and  is  not  dependent  on 
a linear  FM  characteristic.  For  slowly  varying  modulations  and  short  time  delays, 
the  range  response  of  an  IF  replica  correlator  is  given  by 


RL(t)  = R0(T-ir)  + R0(T+Tr)  (2.1) 

where  Rq(t)  is  the  autocorrelation  of  the  transmitted  RF  signal  and  xr  is  the 
reference  delay. 

The  above  result  is  importai  because  it  demonstrates  that  the  range  re- 
sponse of  IF  correlators  can  be  ex...essed  in  terms  of  a statistical  property  of 
the  signal  (autocorrelation)  and  is  therefore  not  dependent  on  the  exact  time  his- 
tory of  the  signal  and  reference.  Following  discussions  will  demonstrate  that 
any  reference  can  be  expressed  in  terms  of  replica  references  and  the  range  re- 
sponse expressed  in  terms  of  autocorrelation  functions;  therefore,  any  IF  correlator 
can  be  analyzed  by  the  statistical  techniques  appropriate  for  IF  replica  correlators. 

For  many  systems  (harmonic  N systems,  for  example)  the  reference  is  derived 
directly  from  the  modulation  circuitry.  Simple  systems  of  this  type  are  usually 
dependent  on  a linear  modulator  characteristic.  If  the  modulator  is  linear,  a 
reference  can  be  generated  which  differs  from  a replica  reference  only  by  a phase 
constant  (due  to  the  fact  that  the  modulator  has  no  information  about  oscillator 
center  frequency  u)q).  The  phase  constant  has  no  practical  efr  t on  the  range  re- 
sponse since  it  appears  as  a doppler  phase  angle.  However,  the  range  response  is 
not  strictly  a displaced  autocorrelation  function  since  the  doppler  carrier  will 
not  peak  at  the  same  point  as  the  range  law.  References  which  differ  from  replicas 
only  by  a phase  constant  will  be  called  quasi-replicas;  the  range  responses  which 
result  from  these  references  will  be  called  quasi-autocorrelation  functions. 

3.  STATISTICAL  DESCRIPTION  OF  SIGNALS  AND  REFERENCES 


This  section  describes  the  signal  and  reference  in  terms  of  0 and  discusses 
the  transition  from  time  to  statistical  averaging,  A graphical  technique  for 
determining  R(0)  is  described  including  a discussion  of  the  procedure  to  be  fol- 
lowed when  the  reference  is  not  a single-valued  function  of  0. 

3.1  Signal  Description 

A constant-amplitude  angle-modulated  carrier  can  be  represented  as 

x(t)  * cos  0 i (t)  (3.1) 

t 

where  0j(t)  ■ / wc(A)dA  and  ojc  is  the  instantaneous  frequency.  It  is  convenient 

to  describe  wc  as  the  sum  of  a constant  center  frequency  u>q  and  a valving  component 
<%(t).  ©i  (t)  then  becomes 

t t t 

0 1 (t)  B / wgdA  + / uJm(A)dA  » ajgt  + ^(AjdA  . (3.2) 
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A delayed  return  signal  can  be  written  as 

x(t-t)  = cosCOjU-t)] 

where  t is  the  delay  time.  0i(t-T)  is  given  by 


(3.3) 


Oj (t-r)  = 


0 


t-T  t-l 

/ at()dX  + / u>m(X)dA 


t-T 


= (Oq t - wot  + / ii>m(X)dA 


(3.4) 


The  IF  signal  is  obtained  from  some  process  (envelope  detection,  etc.)  which  is 
equivalent  to  multiplying  x(t)  by  x(t-x)  and  low-pass  filtering  the  double  fre- 
quency terms  to  obtain  the  lower  sideband.  The  IF  signal  is  then  given  by 


s(t)  = cos[0i(t)  - 6 i ( t— t ) ] 


= cos[(i)QT  + / ^(XjdX] 

t-T 


(3.5) 


Expanding  I m(T)dT  in  a Taylor  series  about  t » 0 gives 

t-T 


/ (um(X)dX  - u)m(t)T  - 
t-T 


£>m(t>^ 


(3.6) 


<em(t)T2 


For  many  ranging  systems,  the  phase  error  represented  by  2 


+ •••  is  small 


and  the  signal  phase  is  adequately  represented  by  the  angle  (u>0  + <%(t))T.  The 
signal  can  then  be  written  as 


S (t)  = COS [tOc  ( t ) T J = COS [u>m( t ) T + W()t] 


(3.7) 


The  signal  is  now  seen  to  be  a function  of  the  angle  u>ct  without  regard  to  the 
time  history  of  The  modulation  angle  0 is  defined  as 


0 = 


(3.8) 


where  tq  is  a chobcr.  reference  delay  (to  be  discussed  later). 
The  signal  can  now  be  written  in  terms  of  0 as 


S(0)  = cos [u>0 t + 0 


(3.9) 
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3.2  Transition  From  Tine  Averaging  to  Statistical  Averaging 

IF  correlator  systems  may  be  represented  by  the  partial  block  diagram  of 
Figure  2(a). 


(A)  Hardware  and 
Real  Time 
Diagram 
s(t) 


(B)  Equivalent 
Statistical 
Math  Model  1 


* RL  - s(t)*r(t) 


s(0)  * 

Statistical 

Averager 

<c>  SSSSEi  !*<m> 

Math  Model  2 

s(0)  » 

Statistical 

-Hf 

Averager 

-►  RL  - E{S(0)  • /?  (0 ,i) } 


-►  RL  =»  E{S(0)-R(6)} 


R(0)  - l P(i|6)/?(0,i) 
i 


Figure  2:  IF  Correlator  Partial  Diagrams 


The  signal  and  reference  waveforms  are  multiplied  and  time-averaged  by  the  filter 
to  yield  the  range  response 


RL  * s(t)r(t) 


(3.10) 


Jn  order  to  treat  periodic  systems  by  statistical  methods,  it  is  useful  to  treat 
time  as  a random  variable  and  all  functions  of  time  as  transformations  of  the 
random  variable  t.  A suitable  representation  of  t is  that  of  a uniformly  distri- 
buted random  variable  with  probability  density 


1, 

T 


p(t) 


T , T 
2 < < 2 


otherwise 


(3.11) 
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Thus  we  may  replace  all  time  averages  by  ensemble  averages: 


T/2  

E{s(t)r(t) } = / s(t)r(t)  ^ dt  * s(t)r(t)  . (3.12) 

-T/2 

This  relation  reflects  an  assumption  of  ergodicity  of  the  processes,  i.e.,  the 
time  average  equals  the  ensemble  average. 

Now  consider  a general  periodic  modulation  consisting  of  a number  of  monotonic 
sections  as  shown  in  Figure  3.  We  denote  the  sections  by  index  i and  their  length 
in  time  by  Tf.  Now,  if  we  define  a function  0 to  be 

0 * gi(t)  = g(<%(t))  * a>in(t)T0  (3.13a) 

where  g is  a functional  of  and  tq  is  a chosen  delay,  and  since  i is  also  a func- 
tion of  t, 

i **  g2(t)  , (3.13b) 


<>(t)  - ( t ) i o 


Figure  3:  Monotonic  Modulation  Sections 


then  there  is  a one-to-one  correspondence  between  the  two  sets  ftl  ,ind  ( ,11  such 
that,  using  the  transformation  of  variables  by  equation  (3.13),  p(t)  maps  into 
the  joint  probability  density  p(0,i)  and  vice  versa.  It  is  seen  from  equation  (3.9) 
that  the  signal  S(0)  is  not  a function  of  i.  The  reference  may  also  be  written  in 
terms  of  0,  but  the  reference,  /?(0,i),  is,  in  general,  also  a function  of  i.  Tt  is 
clear  that  we  may  now  average  with  respect  to  any  statistic  in  which  the  argument 


5 

! 


I 


of  the  average  is  expressed. 

E{s(t)r(t) } - E{S(0)#(0,i)}  (3.14) 

where  we  have  indicated  that  the  signal  and  reference  may  be  expressed  in  either 
domain  by  a transformation  of  variables.  Figure  2(b)  illustrates  this  notation. 

Now 


E{S(0)i?(0,i)}  « /S(G)[//?(o,i)p(0,i)di]d0  . (3.15a) 

Since  the  joint  density  p(6,i)  is  that  of  a continuous  variable  0 and  a discrete 
variable  i,  integration  with  respect  to  i gives 


E{S(0)fl(0,i)}  = /S(0)[£/?(0,i)P(i|o)]p(0)d0 
= /S(0)R(0)p(0)d0 

= E{S(0)R(0) } (3.15b) 


where  R(0)  can  be  considered  a composite  reference  described  by 


R(e)  - I /?(6 ,i)P(i|  6)  . (3.16) 

i 

P (i| 0)  is  the  conditional  probability  of  the  discrete  random  event  i. 


We  may  now  equivalently  show  the  variables  as  in  Figure  2(c).  The  description 
of  S and  R in  terms  cf  0 will  be  used  in  the  statistical  methods  described  in  this 
report.  Equation  (3.15)  shows  that  any  expression  for  R(0)  which  is  correct  over 
the  modulation  interval  bounded  by  p (0 ) will  give  the  same  range  response.  How- 
ever, the  mathematical  form  of  the  results  may  appear  quite  different  since  various 
representations  of  R(0)  which  are  equivalent  over  the  modulation  interval  may  be 
quite  different  beyond  that  interval.  This  point  is  illustrated  in  an  example  to 
follow. 

3. 3 Description  of  R(0) 

3.3.1  Determination  of  R(0)  from  r(t)  and  0 (t ) 


As  noted  previously,  the  relationship  between  the  reference  and  the  modulation 
is  determined  explicitly  by  the  hardware  implementation  and  can  be  written  analyt- 
ically if  desired.  However,  a simpler  graphical  technique  can  be  used  to  determine 
R(0)  when  r(t)  and  0(t)  * ootn(t ) tq  are  kn^wn.  The  linear  modulator  characteristic 
is  assumed  so  that  <%  is  linearly  proportional  to  the  modulation  voltage  vm  as 
given  by  0 * DvnTo.  A visual  display  of  the  functional  relationship  R(0)  can  be  ob- 
tained by  connecting  r(t)  to  the  y-axis  and  0 ( t)  to  the  x-axis  of  an  oscilloscope 
and  observing  the  Lissajous  figure  which  results.  Alternately,  a graphical  transfer 


function  approach  can  be  used  where  6(t)  and  r(t)  are  plotted  and  G(t)  used  to 
map  r(t)  into  R(6).  The  approach  is  illustrated  in  Figure  4 with  triangle  modu- 
lation and  a sinusoidal  reference.  r(t)  is  written 


2tt 

r(t)  * cos  y t 


and  R(6)  is  written  most  simply  as 


R(8)  = cos  (0  - j) 


(3.17) 


(3.18) 


where  the  range  of  0 is  obviously  it,  implying  that  Btq  = 1/2,  where  B is  the  peak- 
to-peak  RF  bandwidth  in  Hz.  More  formal  mapping  procedures  can  be  defined  but  are 
not  necessary  for  analyzing  the  systems  described  here. 

It  is  apparent  from  Figure  4 that  R(0>  does  not  depend  on  time  period  T; 
therefore,  this  system  may  be  period-modulated  as  in  an  FM/FM  type  system  without 
changing  R(0)  or  the  range  response  if  r(t)  and  0(t)  are  generated  by  methods  which 
do  not  change  their  functional  relationship  when  modulated.  Multiple-output  func- 
tion generators  which  may  be  frequency-modulated  represent  such  systems. 


o(t)  = u>m(t)  t o 


Figure  4:  Graphical  Mapping  of  r(t)  to  R(0) 


\ 


For  the  previous  example,  a given  value  of  0 corresponded  to  a single  value 
of  r(t).  A replica  reference  always  has  this  one-to-one  relationship;  however, 
an  arbitrary  modulator/deoodulator  system  may  not  produce  such  a unique  (i.e., 
one-to-one)  relationship.  For  such  cases,  we  perform  the  mapping  procedure  for 
the  separate  monotonic  sections  of  the  modulation  (a)m(t)  or  0(t))  to  produce  a set 
of  references  {/?(0, i)}  where  i refers  to  the  modulation  section.  To  illustrate 
this,  consider  the  previous  example  with  a reference  phase  of  90  degrees  (Fig- 
ure 5).  The  range  response  of  this  system  will  be  zero  since  the  reference  is 
not  phased  properly  (time  domain)  with  respect  to  the  modulation.  It  is  seen 
that  the  first  section  of  °(t)  maps  positive  values  of  r(t)  into  R(0,1)  while 
the  second  section  maps  negative  values  of  r(t)  into  R(8,2).  To  carry  out  the 
analysis  indicated  for  cases  of  this  type,  a single,  composite  R(0)  is  determined; 
the  functions  R(0,i)  are  simply  averaged  with  respect  to  i to  yield  R(0).  Ac- 
cording to  Section  2.2,  we  form  this  average  as 


R(0)  - l P(i|0)fl(0,i) 
i 


(3,19) 


<D  CD 


C4 

* 

CD 

* 

as 


e(t)  * W|n(t)To 


Figure  5:  The  Multi-Valued  Mapping  Case 
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It  is  easy  to  show  that,  in  general,  the  (discrete,  <’onditional)  probability  of 
the  ith  modulation  section  is  found  as 


I 


& 


K 


5' 


IT 


1 


! 

5 


s 

t 

v 

& 


P (i | 0 ) 


lSj-l(e)| 

Z I Sj“ 1 (Q) I 
j 


(3.20) 


where  Sr  * -£r  0(t)  is  the  slope  of  the  kth  section  of  0(t)  « wm(t)io  at  the  ordi- 
nate 0 (and  is,  in  general,  a function  of  0);  the  summation  of  the  denominator  is 
taken  over  those  sections  which  pass  through  the  ordinate  0. 


For  the  current  example  IsJ  = |s2|  and  thus  P (X | 0)  * P (2 | 0)  = 1/2.  As  ex- 
pected 

i 

R(0)  = - /»(0,1)  + j /?(0,2)  - 0 . * (3.21) 

Since  the  composite  reference  is  zero,  the  system  will  have  no  range  response  as 
expected. 

3.3.2  Choice  of  tq 


The  constant  t0  used  to  define  0 *=  (%t,  is  a scaling  parameter  which  deter- 
mines the  periodic  interval  for  expanding  R(0)  in  a Fourier  series,  tq  is  chosen 
so  that  R(0)  may  be  represented  in  the  most  concise  cosine  series  that  represents 
R(0)  in  the  modulation  interval.  This  is  generally  accomplished  by  defining  the 
smallest  angular  interval  over  which  the  reference  is  periodic  to  be  from  -tt  to  v 
radians  (Figure  6) 


"(t)  * ^(t) 1 n 


Figure  6:  Illustration  of  Assignment  of  Tq 


Note  that  the  total  modulation  angle  2tiBtq  may  be  less  than  or  greater  than  2:;. 
In  some  cases  it  may  be  helpful  and  completely  valid  to  extend  the  definition  of 
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R(0)  beyond  the  modulation  interval  to  obtain  a simple  representation  of  R(6). 

The  ratio  K of  the  interval  chosen  as  ±n  to  the  modulation  interval  ±tiBtq  defines 


to  as 


JL_ 

T°  " KB 

A replica  reference  for  time  delay  tr  ia  defined  by 


(3.22) 


i 0Tr 

R(0)  = S(0)|  = cos(o)0tr+  — ) . (3.23) 

T“Tr  TQ 

It  is  evident  that  for  this  case  tq  should  be  chosen  as  tq  ■ ir.  A poor  choice 
of  tq  will  still  yield  the  correct  range  response  but  will  require  that  more  terms 
be  sunned  to  give  the  same  results.  The  examples  to  follow  will  clarify  the  choice 
of  t0. 

3.3.3  Series  Expansion  of  R(0) 

To  facilitate  the  following  analyses  F(0)  is  expanded  in  a harmonic  cosine 
series.  Having  defined  the  periodic  interval  as  described  in  the  previous  section, 
R(0)  can  be  written  as 

! 

R(0)  - l Cn  co8(n0  + <|>n)  ; -itBto  < 0 < ttBtq  (3.24) 

! n-0 

j 

j where  B is  the  peak-to-peak  frequency  deviation  in  Hz, 

| Cn  - K2  + Bn2  * ♦n  ■ tan"1  (~^) 

1 1 71  l 11 

J A0  - / R(9)d0  , An  - - / R(0)cos  n0d0 

j ”TT  “IT 

1 ^ 

J Bn  * — / R(0)sin  n0d0 


The  series  is  the  sum  of  quasi-replicas. 

A series  of  harmonic  replicas  would  be  written  as 

00 

R(0)  * l Cn  cos[nO  + nojoio]  • (3.25) 

n-0 

The  quasi-replica  is  distinguished  from  the  replica  in  that,  in  general,  4>n  / n(.>oT0 
since  quasi-replicas  are  independent  of  center  frequency,  wq. 
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4 . STATISTICAL  ANALYSIS  OF  IF  CORRELATORS 

Equation  (3.15)  provides  the  basis  of  statistical  analysis.  This  statistical 
average  may  be  expressed  variously  as 


RL(t ) - E( S(0)/7(0,i)  ) 


- /S(0)[}>(0,i)P(i| 0)]p(0)d0 

i 

= /S(<i)R(O)p(0)dO 

- ^P(i)[/S(0)/?(0,i)p(0 1 i)do] 
i 

- [P(i)RLi(T> 
i 


(4.1a) 


(4.1b) 


The  average  expressed  as  in  Equation  (4.1a)  implies  that  a composite  reference 
R(0)  is  computed  as  the  statistical  average  (over  i)  of  the  reference  H(Ut i)  cor- 
responding to  the  modulation  sections.  This  reference  is  then  used  in  computing 
the  first-order  average  with  respect  to  6.  A secoi  d expression  for  the  average, 
equation  (4.1b),  computes  the  contribution  of  the  ith  modulation  section  (with  the 
corresponding  reference)  and  then  averages  the  resulting  component  range  responses 
RLi(t).  In  this  report  we  choose  the  former  approach  because  the  range  response 
may  then  be  expressed  in  terms  of  the  signal  autocorrelation  function. 

The  form  of  Equation  (4.1a)  shows  that  systems  are  equivalent  (i.e.  yield 
the  same  range  response)  as  long  as  the  product  S(G) ’R(0) *p(0)  remains  the  same. 
This  means  that  any  combination  of  FM  weighting  (reflected  in  p(Q)),  AM  weighting 
of  the  signal  envelope  (reflected  in  S(0)),  or  weighting  in  the  "receiver"  (re- 
flected in  R(6))  may  be  used  to  achieve  a desired  range  response. 

For  the  constant-amplitude  FM  signal  case  with  R(0)  expanded  as  a cosine 
series  the  (normalized)  range  response  is 


RL(t)  « 2/cos(0  — + u)q t ) J Cn  cos(n0  + <j>n)p(6)d0 
T°  n 

* T Cn{  / cos[G(—  + n)+(u>oT+<|>n)]p(0)d0+/  cos[0(— — 
n T°  T° 

= l Cn{cos(u>0T  + <pn)  / cos[0(~*  + n)]p(6)dG 
n 0 

- sin(o)0T  + <j>n)  f sin[o(-^-  + n)]p(0)dO 

+ cos (uiqt  - 4>n)  / cos[0(— — n)]p(0)d0 

T0 

- sin(o)QT  - tn)  / sin[0(— n)]p(0)d0} 

T0 


(4 . 2n) 

n)+(u)0T-<{>n)]p(»)do ) 

(4.2b) 


(4.2c) 
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Since  sin[8( — ± n)]  is  an  odd  function  of  0,  if  the  power  spectrum  is  even 
T0 

(i.e.  p(0)  is  even),  the  second  and  fourth  integrals  of  Equation  (A. 2c)  are  zero, 
simplifying  the  expression  to 

RL(t)  » £ Cn{cos(u)oT  + <j>n)  /cos[o(— + n)]p(0)d6 
n 0 

+ cos(o)0T  - (|>n)  / cosl«(~-  - n)]p(0)dO)  . (A. 3) 

To 


The  integrals  are  recognized  as  Fourier  cosine  transforms  of  p(0): 

/ cos[0 (~  ± n)]p(0)d0  = p0(i  ± nx0)  (A. A) 

To 

where 

Po(t)  " F 1(F0 (<*)))  , 

Fo(u)  ■ 2tttop(8)  with  the  substitution  0 * ojtq 

and  Fq  is  the  signal  power  spectrum  translated  to  baseband.  Thus  the  range  re- 
sponse becomes 

RL  = E{SR}  * l Cn[po(T-nT0)  cos  (u)0T-*<f>n)  + p0(t+nT0)  cos  (w0T+<j>n)]  • (A. 5) 

n 

Equation  (A. 5)  gives  the  range  response  in  terms  of  pp,  doppler  phasing  <f>n,  and 
weighting  coefficients  Cn.  Po  is  related  to  the  envelope  of  the  power  spectrum 
of  the  transmitted  signal  (translated  to  baseband)  through  the  Fourier  transform; 
{Cn}  and  Hn}  are  derived  from  R expressed  as  a harmonic  series  in  0. 

It  is  important  to  realize  that,  as  long  as  the  relationship  of  R to  0 re- 
mains unchanged  and  the  RF  power  spectrum  remains  the  same,  then  the  range  re- 
sponse does  not  change  (note  the  parameters  of  Equation  (A. 5)).  Therefore,  the 
statistical  analysis  technique  may  be  used  with  equal  facility  for  both  harmonic, 
random,  and  FM/FM  type  systems. 

5.  EXAMPLES  OF  MAPPING  AND  ANALYSIS 


The  examples  of  this  section  are  designed  to  interrelate  various  concepts, 
to  give  insight  into  the  statistical  analysis  approach,  to  contrast  that  approach 
with  time  analysis,  and  to  illustrate  the  basic  method.  The  examples  will  follow 
the  following  analysis  procedure. 

Step  1:  Determine  R(  0) 

(A)  By  knowledge  of  r(t)  ■ 0(t)  * <%(t)io,  determine  R(0)  or  [H('\ 0). 

(The  graphical  mapping  procedure  will  be  used  in  the  examples.) 

(B)  Average  U?(0,i)}  as  shown  in  Section  3 (if  the  mapping  is  not  one- 
to-one)  to  yield  R(8). 
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Step  2:  Expand  R(fl)  into  a harmonic  cosine  series. 

Define  tg  to  yield  tne  most  concise  series. 


R(0)  s 1 C„  cos  (n0  + 4>n) 
n*0 

Step  3:  Write  the  range  response. 

Assuming  at*  even  power  spectrum,  the  response  will  be 

oo 

RL  = l Cn[p0(T-nT0)  cos  (u>0T  - 4in)  + p0(T+nT0)  cos  ('.igt  + d>n) 3 
n=0 


(5.1) 


(5.2) 


The  mapping  of  r(t)  into  R(0)  is  central  to  the  analysis.  The  following  examples 
illustrate  this  mapping  and  will  clarify  other  features  of  the  statistical  approach 
to  the  analysis.  In  most  examples,  the  time  histories  of  r(t)  and  {vm,  <%,  wmXg-O} 
are  depicted  showing  the  mapping  into  R(Q, i)  and  thus  RvO).  The  probability  den- 
sity p(G)  is  also  illustrated,  since  it  determines  the  envelope  of  the  transmitted 
power  spectrum  and  thus  Po(t). 

Example  1:  IF  Replica  Correlators 

The  replica  correlator  system  generates  a reference  which  is  an  exact  copy 
of  the  IF  signal  for  t - tr,  where  xr  is  the  delay  for  which  the  response  is  de- 
signed to  peak.  Thus 

r (t)  * s(t)  I = cos  (oc(t)Tr  « cos  [wm(t)Tr  + w0xr] 
i r 


cos  [o (t)  ~ + U)0Tr] 


(5.3) 


We  now  follow  the  steps  for  finding  the  range  response. 


Step  Is  Write  R(0). 


R(0)  = cos  (0  — + 0)gTr) 


Step  2:  Expand  R(0)  in  a series. 

As  mentioned  previously,  our  choice  of  tq  is  flexible,  and  is  chosen  to 
give  the  most  concise  cosine  series  for  R(8).  We  therefore  choose  ig  = Tr  to 
give 

oo 

5(0)  = l Cn  cos  (n0  + <t>n) 
n=0 
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Figure  7:  Hard-Limited  Replica  Reference 


Step  3:  Using  (5.2)  the  range  response  is 


RL  = l Cn[p0(T-ntr)  cos  wo(T-nTr)  + p0(x+nTr)  cos  uin(x+mr)] 

= l Cn[R0(T-nrr)  + R0(T+nTr)] 

It  is  important  to  note  that  each  component  of  RL  has  an  identical  envelope 
which  depends  only  on  the  envelope  of  the  transmitted  power  spectrum;  that  is, 
it  depends  on  the  probability  density  of  the  modulation  and  not  specifically 
on  the  modulation  waveshape. 

We  see  that  the  range  response  of  a system  which  employs  a replica  which 
has  been  distorted  by  any  general  zero-memory  nonlinearity  will  consisc  of  the 
same  translated  autocorrelation  responses  as  in  the  above  expression;  only  'Cn' 
will  vary. 

Example  3:  Quasi-replica  System  Using  Sinusoidal  Modulation 

This  example  is  interesting  because  it  demonstrates  that,  while  the  time- 
domain  or  spectral  representation  of  the  reference  is  complex,  its  functional 
description  is  quite  simple.  The  reference  R(0)  is  assumed  to  be  identical  to 
a replica  except  for  phase.  Also,  the  reference  will  be  assumed  to  be  derived 
from  the  modulation  voltage  rather  than  the  RF  carrier  (i.e.,  there  are  always 
4 reference  cycles  per  modulation  cycle  regardless  of  RF  bandwidth).  For  this 
reason,  the  reference  will  match  the  returned  signal  at  a specific  value  of  B 
rather  than  x alone  as  would  be  the  case  for  a replica  derived  from  the  RF 
carrier.  This  distinction  holds  for  all  systems  which  derive  the  reference  from 
the  modulation  voltage  rather  than  the  RF  carrier. 
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Figure  8!  Quasi-Replica  Reference  for  Sinusoidal  Modulation 


Step  1:  The  non-uniform  time  reference  r(t),  when  mapped  through  0(t),  yields 

R(0)  which  is  purely  sinusoidal. 

Step  2:  Choosing  the  smallest  periodic  interval  defines  Bt0  *»  2 or  x0  » 2/B. 
R(0)  is  described  as  R(0)  - con  0.  Using  equation  (5.1)  R(0)  is  now  written  in 
its  simplest  series  with  Cj  * 1,  Cn  * 0 for  n i 1,  and  <|>n  « 0. 

Step  3:  Using  equation  (5.2),  the  response  is  written 

RL  “ [pu(t_to)  + PO  (t+to)3  cos  k>0T 

where  the  autocorrelation  envelope  is  found  from  the  power  spectrum  to  be 


Po(T)  * 

Since  tq  ■ 2/B,  the  range  response  becomes 


RL  - [J0(it[Bt-2])  + J0(*[Bt+2])]  COS  UIqT 
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A time  analysis  of  the  product  s(t)r(t)  would  have  been  much  more  difficult 
since  the  contributions  due  to  all  harmonics  must  be  summed  to  give  the  range 
response.  We  see  in  this  example  the  simplicity  of  statistical  analysis — the 
transformation  to  the  0-domain  transforms  all  replicas  or  quasi-replicas  to 
simple  sinusoidal  form  no  matter  how  formidable  they  appear  in  time  representa- 
tion. 


ri 

■a 


Example  A:  Triangular  Modulation  with  Sinusoidal  Reference 

Figure  9 illustrates  the  system  waveforms.  This  system  may  be  called  an 
N*2  harmonic  system  3ince 

r (t)  = cos  [2  • Oy-)t] 


Steps  1 and  2:  Since  0(t)  is  linear  and  the  reference  r(t)  is  sinusoidal,  R(o) 

is  also  sinusoidal.  The  smallest  periodic  interval  is  the  total  modulation  in- 
terval so  that  Sto  = 1 or  tq  = 1/B.  Then 

R(0)  = -cos  0 

C i — —1 , Crj  = 0 * n ^ 1 


*n  = 0 


n 


W 


ir>>»fwmgsraa  wwa.*  waw^jca* . »,,  w,  ^ 


Step  3:  Using  equation  (5.2),  the  range  response  becomes 

RL  ■ Ci [p0 (t-t0)  + Po(t+t0)3  cos  u0t 


Since  the  power  spectrum  is  uniform, 


and 


Po(T) 


sin  ttBt 
nBx 


I! 


■§ 


sin  tt(3t-Btq) 

RL  * -[ C08  U)qT 


it(Bt-Btq) 
sin  7t(Bt+Bt0) 

71  (Bt+Btq) 


cos  ojqt] 


Substituting  tq  * 1/B  yields 


„„  rsin  tt(Bt-I)  . sin  ti idti-i;  n 

RL  ' "I' ' TTuV^T"  008  “ot  + ~(gl+iT~  cos  “°t] 


sin  tt(Bt-H) 


I 


This  result  can  be  verified  by  time  analysis. 

Example  5:  Triangular  Modulation  with  Squarewave  Reference 


1 


Steps  1 and  2:  This  example  corresponds  to  Example  4 with  the  reference  defined 

as 

r(t)  = sgn[cos  ~ t]  . 

Thus  R (0 ) * sgn  [-cos  0],  and  xg  = 1/B  as  in  Example  4.  Using  equation  (5.1), 
the  squarewave  in  0 is  expanded  as 


R(e)  = l Cn  cos  n0 
n=0 


r 4 , nir 

Cn  = sin  -=- 

11  n7t  2 


Step  3:  Using  equation  (5.2)  the  range  response  is  written  directly  ns 


CO 

v ^ rsin  ii(BT-n)  . sin  n(Bx+n)i 
RL  - cos  u„t  l C„  [-T(-bV:-T  l * nTBr+n) 


n=l 


These  results  are  the  same  as  those  of  Example  4 except  for  the  addition  of 
range  responses  with  displacements  and  amplitudes  according  to  the  terms  of  the 
harmonic  series  for  R(0). 
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Example  6:  Elimination  of  Third  Harmonic  by  Reference  Phasing 


Squarewave  references  are  more  economical  to  generate  than  sinusoidal 
references,  but  they  add  undesirable  range  responses  due  to  their  harmonic  con- 
tent. Generally,  harmonics  higher  than  the  third  are  not  troublesome  since  sig- 
nal losses  combined  with  lower  amplitudes  reduce  these  responses  to  a tolerable 
level.  Filtering  the  IF  signal  to  eliminate  components  corresponding  to  ref- 
erence harmonics  has  been  used  [5],  but  this  technique  requires  more  circuitry 
chan  the  phase-shift  technique  of  this  example.  Shifting  the  reference  with 
respect  to  the  modulation  by  30  degrees  achieves  a third-harmonic  reference 
phase  of  90  degrees,  which  gives  zero  third-harmonic  range  response.  All  other 
responses  present  will  be  attenuated  by  0.866. 


csi 

* 

<T’ 


CD 

w 

a 


2 11 B 1 1 
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Steps  1 and  2:  As  shown  in  Figure  10,  the  reference  averages  to  zero  in  the 

double-valued  region  (since  the  weightings  of  /?(0,1)  and  i?(0,2)  are  both  1/2). 
The  reference  may  be  expanded,  with  Btq  * 1,  as 

OO 

R(0)  = £ Cn  cos  n0 

n=0 


C 


n 


cos 


nTT 

6 


n odd 


0 


n even 


Note  that  C3  **  0,  eliminating  the  third  harmonic  of  the  composite  reference,  and 
thus  the  third  harmonic  response. 

Example  7;  Ji  Ranging  System 

This  example  demonstrates  the  corroboration  of  functional  analysis  with 
conventional  analyses  techniques  for  the  system  known  to  give  the  range  response [9] 

RL  = -2Ji(ttBt)  sin  ujqt 

This  example  also  demonstrates  that  the  range  response  can  be  expressed  in  various 
equivalent  forms  since  any  expression  for  R(0)  which  is  correct  in  the  modulation 
interval  bounded  by  p (9 ) can  be  used.  Figure  11  shows  the  sinusoidal  modulation 
and  demodulation  in  the  time  domain,  along  with  the  linear  function  R(0)  and  p(G) 
(which  describes  the  envelope  of  the  transmitted  power  spectrum) . 


* 


2ttBt  0 

I 


Figure  11:  Ji  Ranging  System 


Step  1:  R(0)  is  linear,  and  with  axes  labeled  as  in  Figure  11. 

R(0)  * - 6 

71 

Step  2:  It  is  convenient  to  let  the  range  of  the  angle  modulation  by  ± tt/2  (or 

Btq  ■ 1/2)  which  implies  half-wave  symmetry.  R(0)  is  then  a triangular  wave  which 
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has  been  defined  over  the  interval  ± it  (beyond  the  limits  of  the  modulation)  to 
simplify  the  series  representation.  We  then  describe  R(G)  in  a harmonic  series. 

00 

R(0)  = l Cn  cos  (n0  - y) 
n*l 


where 


8 . nir  . , it 

— sin  -y  and  *n  ■ “ 2 


Step  3:  The  range  response  is 


RL  - l Cn[po(T-nTo)  cos  (uiqt  + y)  + po(T+nio)  cos  (u)0t  - y)3  . 


For  a sinusoidally  modulated  carrier  P0(t)  is  given  by 


p0(t)  - Jq(itBt)  , 


and  since  Bt0  ■ 1/2, 


RL  - l Cn{Jo[^(BT  - y)3  COS  (w0T  + y) 
n-0  1 1 

+ J0[tt(Bt  + y)3  cos  (o)0T  - "■)} 

- -sin  U)0T  l Cn{ J0 [ tt (Bt  - y)3  - J0[*(Bt  + |)]}  . 

It  is  interesting  to  see  the  range  response  in  terms  of  displaced  Jq  responses. 
To  be  confident  that  these  results  are  correct  we  may  check  to  see  if  J)  is  in- 
deed expressed  in  terms  of  the  Jq  series.  Taking  a few  as  three  terms,  and 
checking  point-by-point,  we  find  that,  to  a good  approximation, 

2J j (ttBt ) : 4r  U0I>(Bt  - y)3  - J0I>(Bt  + y)] 

- £ J0[tt(Bt  - |)]  +y  J0r»(BT  + f)3 

+ -k  Jo["(Bt  - y)3  - J0I><Bt  + y)3)  . 


The  series  is  approximately  equal  to  2Ji(nBt)  for  Bt  < 10.  Presumably,  if  the 
series  is  not  truncated,  the  sum  is  identical  to  2Ji(ttBt).  It  is  seen  that  many 
(perhaps  useful)  series  expansions,  such  as  this  one  relating  Jq  and  J),  may  be 
generated  using  such  analyses. 
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It  was  demonstrated  earlier  that  any  series  expansion  which  describes  R(6) 
in  the  modulation  interval  will  give  the  same  range  response.  This  may  be  veri- 
fied for  this  example  by  defining  Bt0  = 1 which  then  describes  R(6)  as  a sawtooth 
wave  whereas  defining  Btq  * 1/2  described  R(6)  as  a triangular  wave.  The  series 
for  a sawtooth  converges  less  rapidly  than  the  triangle  series  and  consequently, 
more  terms  must  be  summed  to  obtain  the  range  response. 

A much  simpler  description  is  found  by  considering  R(0)  to  be  a small  por- 
tion of  a large  sinusoid  with  amplitude  A as  given  by 


R(0)  = A sin  0 

Assuming  sin  8 « 8,  A is  found  by  setting  R(0)  * 1 when  0 * ttBtq.  R(0)  then 
becomes 


R(9)  -fia-i.  tm  . 

ttBtq  ttdTq 

Substituting  Ci  = and  $i  = -v/2  gives  the  range  response  as 

sin  cjqt 

RL  - — Lp0(t+t0)  ■ po(t-t0)J 

The  choice  of  t0  can  be  made  arbitrarily  small  (consistent  with  estimating 
sin  0 » 0)  and  po(t-tq)  - Po(T+to)  can  written  as  a differential  as 


lim  [pc(t+t0)  - Po(t-tq)]  * 2t0Po’(t) 
t0-K) 

The  range  response  is  then  given  by 


po’Ct) 

RL  = 2 sin  <dqt  — — — 

This  result  holds  for  all  IF  correlators  which  use  the  modulation  as  the  IF 
reference. 

For  this  particular  case  pq(t)  = Jo('nBt).  Differentiating  with  respect  to  t 

gives 

P°’(t)  = [ J0  (ttBt  ) ] * -ttBJ  j (uBt  ) 

The  range  response  is  then  equal  to  the  result  known  to  be 

RL  =~2J i (ttBt)  sin  wqT 

For  this  particular  case,  the  equivalence  of  different  expressions  for  the 
range  response  could  be  shown  to  be  identical;  however,  in  many  cases,  this 
equivalence  may  not  be  easily  demonstrated. 
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Example  8:  Hybrid  Modulation  Systems 


The  systems  described  in  this  example  are  designed  to  fulfill  three  objec- 
tives: 

(1)  Provide  a range  response  with  good  definition  by  controlling  the 
envelope  of  the  trenemltted  power  spectrum; 

(2)  Use  quasi-replica  correlation  techniques  to  obtain  displaced 
responses  without  using  a delay  line  to  produce  the  reference; 

(3)  Synthesize  a system  which  may  be  period-modulated  (FM/FM)  to 
produce  a wideband  IF  in  order  to  reduce  the  frequency  ambi- 
guities inherent  in  harmonic  systems. 

Figure  12  illustrates  two  different  hybrid  modulation/demodulation  techniques 
which  produce  the  same  power  spectra  and  range  responses.  Both  techniques  are 


9(t)  * Wm(t)T0 


r(t) 


r(t) 

:i;AAAAAAA^ 


Figure  12:  Hybrid  Modulation  Systems 


* ^F"A?Wd'»^^^^a!'»Ji!lf>JSW  •'  'JWJBWAJWLA1.?, 


implemented  using  sawtooth  or  triangle  waveform  generators  which  drive  binary 
divider  chains  to  produce  a hybrid  modulation  consisting  of  the  sum  of  analog 
and  digital  variables. 

Steps  1 and  2:  As  seen  in  the  figures,  when  Btq  » 7/2  the  reference  is  written 
as 


R(0)  * cos  0 

Step  3:  Thus  the  range  response  becomes 

RL  = [po(*-t<$)  + p0(T+To)]  c°s  “O7 


Observing  that  the  envelope  of  the  power  spectrum  is  the  superpositioning  of 

component  rectangles,  so  that  the  Fourier  transform  of  the  power  spectrum  is  the 

sum  of  the  envelopesof  the  Fourier  transforms  of  each  component  rectangle,  we  may 

write  c „ _ _ 

. /5ttBtx  . ,3ttBtx 

( ’ ,sl"(  ’ ’ 


16 


(ttBt) 


7 


+ 1 


(3»|T) 


. ,iBTx 
. sin  (-=-) 

+ — — 

7 /TtBt, 


Substituting  Bt0  - 7/2  gives 


RL  - — cosuqT 


78inn(BT-  y)  5siny^Bt-  |-)  3siny-(Bi-  |-)  lsiny(BT-  y) 
+ — + — + 


it(Bt-  -j) 


7. 

T(B,-  y) 


7(Bt*  T> 


7alnn (Br+-  y)  5ain“~(BT+  y)  3sirry-(BT+  y)  lsiny(BT+  y) 
••  "4“  n . t 4" 


tt(Bt+  j) 


^(Br+  |)  -y(Br+  |) 


y(Bi+  {) 


Of  course,  the  use  of  a squarewave  reference  would  result  in  harmonic  range  re- 
sponses according  to  the  harmonics  of  the  squarewave  as  in  Example  5. 

Example  9:  Nulling  of  Sidelobe  Responses  in  Hybrid  Systems 

Example  5 illustrates  the  problem  of  sidelobes  in  the  range  response  of 
squarewave  demodulated  systems  with  uniform  modulation  probability  density.  These 
sidelobe  responses  are  due  to  (1)  the  natural  sidelobes  of  a sin  x/x  shaped  re- 
sponse (caused  by  the  uniform  power  spectrum),  (2)  the  additional  harmonic  re- 
sponses due  to  squarewave  references  R(0),  and  (3)  the  overlap  of  the  positively 
and  negatively  displaced  response  shapes  which  is  significant  when  nBf  is  small. 
Figure  13  illustrates  a simple  system  for  which  the  modulation  is  the  sum  of 
square  and  triangular  waves  and  the  reference  is  the  squarewave.  Ttu  se  waveforms 
can  be  generated  using  economical  integrated  circuit  function  generators  which 
can  be  frequency  modulated  to  obtain  FM/FM  type  systems.  The  ratio  of  square  to 
triangular  wave  in  the  modulation  may  be  adjusted  to  place  a null  in  the  most 
undesirable  sidelobe  as  described  in  reference  [6]. 


Figure  13:  Hybrid  System  With  Null  at  r = 1//B 


In  the  figure  we  have  shown  pi  and  p2,  corresponding  to  the  component  spectra, 
along  with  the  respective  component  references,  /?( 0,1)  and  /?( 0,2).  The  analysis 
may  be  completed  in  two  ways: 

(1)  The  composite  reference  R(0)  and  the  composite  spectrum  corresponding 
to  p (0)  may  be  used  to  obtain  an  expression  for  RL  terms  of  p0,  {Cn } , and  {<|>n}. 

(2)  The  range  response  may  be  regarded  as  the  superpositioning  of  the  re- 
sponses RLi  due  to  /? (0 , 1 ) and  p(o|l),  and  RL2  due  to  /?( 0,2)  and  p(o|2).  Results 
will  not  be  of  the  same  form  as  in  (1),  but,  rather,  in  terms  of  sin  x/x  re- 
sponses with  doppler  carrier  frequencies  of  usq  + ttAB  and  u>0  - nAB.  These  com- 
ponent responses  are  written 

RLi  * pg  (t)  cos  (tog  + ttAB)t 

RL2  = Pg  (t)  COS  (ug  - ttAB)t 

and 

rn(T)  = 

' it(B-AB)t 

Thus  the  total  response  is  written 
RL  = RLj  + RL2 

* Pg(T>[cos(wg  + ttAB)i  - cos (uig  - tiAB)t] 

* Pg(x)  2 sin  wgT  [sin  ttABt] 


The  term  in  brackets  produces  nulls  at  a t = m/AB  (m  * 0,  ±1,  ±2,  •••);  there- 
fore, we  adjust  AB  so  that  the  null  for  m ■ 1 occurs  at  the  point  tnujj  of  the 
most  undesirable  sidelobe  peak: 


Tnull 

This  form  of  RL  shows  clearly  the  nulling  term. 

This  example  illustrates  that  expressing  the  range  response  in  terms  of  p0(x) 
by  determining  the  composite  R(0)  is  not  necessarily  simpler  than  summing  the  range 
responses  from  the  individual  components.  In  this  particular  case,  the  expression 
for  pq  may  be  of  a form  which  provides  little  insight;  certainly  R(0)  is  much 
harder  to  express  than  #(0,1)  and  #(0,2). 

6.  HYBRID  MODULATION  SYSTEMS  AND  TEST  RESULTS 


The  following  systems  illustrate  synthesis  of  desirable  responses  using  the 
analysis  techniques  of  this  report  (figures  14-30) . Note  that  most  of  these 
hybrid  systems  are  constructed  so  that  they  are  easily  period-modulated  to  achieve 
FM/FM  systems. 

Included  in  each  test  example  will  be  a block  diagram  to  show  the  implemen- 
tation, the  graphical  analysis  diagram  showing  modulation  and  demodulation  time 
waveforms,  as  well  as  a functional  description  of  the  reference  and  the  power 
spectrum.  The  resulting  range  responses  obtained  by  measurement  and  by  computa- 
tion are  also  shown.  Measurements  were  made  employing  a variable-length  trans- 
mission line  which  introduces  an  approximate  1/t  attenuation  factor  in  the  results. 
The  line  loss  characteristic  is  shown  in  Figure  30. 

The  first  system  is  a very  simple  system  which  may  be  used  for  wideband  IF 
systems.  The  response  consists  of  superpositioned  sin  x/x  responses  at  amplitudes 
and  displacements  given  by  Equation  (5.2)  using  square-wave  parameters.  The 
second  system  uses  the  same  simple  block  diagram  as  System  1,  but  the  components 
of  modulation  are  added  in  different  amounts,  changing  both  the  functional  rela- 
tionship R(0)  and  the  power  spectrum.  The  effect  is  to  introduce  a null  at  ap- 
proximately Bt  = 5/2,  the  peak  of  the  highest  sidelobe  of  System  1.  System  2 is 
an  implementation  of  the  system  described  in  Example  9.  Systems  3 and  4 corre- 
spond to  the  two  hybrid  systems  of  Example  8.  System  4 can  be  easily  period  modu- 
lated by  frequency  modulating  the  function  generator.  System  3 could  have  been 
implemented  differently  in  order  to  allow  period  modulation.  The  3rd  harmonic 
response  of  Systems  3 and  4 can  be  nulled  by  the  technique  used  for  System  2. 

The  range  response  depends  upon  the  power  spectrum  and  the  functional  rela- 
tionship of  the  modulation  and  demodulation;  thus,  two  principles  are  involved 
in  the  synthesis  of  systems  such  as  those  of  this  section: 

(1)  The  envelope  of  the  power  spectrum  has  the  shape  of  the  pr^  ' «bllity 
density  of  the  modulation  p (0 ) . This  density  is  determined  directly  by  the 
probability  density  of  the  modulating  voltage  p(vm)  if  the  modulation  charac- 
teristic is  linear.  We  may  design  the  spectrum  (and  thus  pq)  by  designing  the 
modulation.  Note  that  there  are  many  modulations  which  will  yield  a desired 
p(0).  (See  Example  8.) 
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(2)  The  functional  relation  R(")  may  be  determined  as  in  Section  4,  given 
the  modulation  and  demodulation  time  functions.  However,  in  synthesis,  we  wish 
to  specify  the  relation  R(«)  to  produce  a certain  response.  The  relation  of  R 
to  o * Dt  ()Vm  may  be  made  inherent  in  hardware  which  operates  on  vm  to  give 
r(t)  * g[  Vm(t)j.  Alternately,  we  may  map  from  the  desired  R(0)  through  u(t)  to 
yeild  an  r(t)  which  is  always  correct,  but  which  may  not  be  conveniently  generated. 
If  r(t)  is  not  conveniently  generated,  alternative  modulations  0 ( t ) with  the  de- 
sired p.d.f.  can  be  used  to  attempt  simplication  of  r(t). 

It  is  clear  from  previous  examples  that  there  exists  freedom  in  the  choice 
of  o(t)  and  r(t)  to  yield  a desired  p(o)  and  R(u).  In  the  systems  described 
here,  a well-defined,  displaced  range  peak  is  desired.  We  therefore  choose  a 
modulation  0(t)  to  approximate  a desired  spectrum  shape  corresponding  to  p(o)  and 
a reference  r(t)  which  maps  approximately  into  a quasi-replica  for  R(e).  Since 
p(0)  and  R(0)  do  not  uniquely  specify  0(t)  and  r(t),  the  designer  must  use  in- 
sight and  a knowledge  of  the  limitations  and  economics  of  available  hardware  to 
determine  an  intelligent  design. 

7.  CONCLUSIONS 


A statistical  analysis  technique  for  IF  correlator  systems  has  been  developed 
under  the  assumptions  of  slowly  varying,  high-index  frequency  modulation  and  small 
propagation  delay.  The  technique  involves  a transformation  of  variables  which  al- 
lows expression  of  the  reference  as  a function  of  the  modulation.  It  is  shown 
that  the  range  response  may  be  expressed  in  terms  of  the  functional  relationship 
of  demodulation  to  modulation  and  the  envelope  of  transmitted  spectrum,  which 
(under  the  high-index  FM  assumption)  is  proportional  to  the  first-order  prob- 
ability density  of  the  modulation.  Because  the  analysis  is  based  solely  upon  modu- 
lation statistics  and  the  functional  relation  R(0),  randomly  modulated  FM/FM  type 
systems  may  be  easily  analyzed  (providing  the  functional  relation  and  modulation 
statistics  may  be  determined),  whereas  time  analysis  of  such  systems  is  often  im- 
practical. The  approach  is  valid  for  all  classes  of  IF  correlators  including 
harmonic  systems  and  replica  correlators  and  unifies  the  treatment  of  IF  systems. 
Systems  which  may  appear  to  be  quite  different  in  the  time  domain  may  indeed  be 
shown  to  be  equivalent  using  statistical/functional  analysis.  The  technique  yields 
the  range  response  in  terms  of  a sum  of  weighted  quasi-autocorrelation  functions 
which  have  been  translated  along  the  t axis.  These  are  determined  solely  bv  the' 
probability  density  of  the  modulation.  The  displacements  and  weightings  depend 
on  the  coefficients  of  a harmonic-series  expansion  of  the  function  R(0).  The 
statistical  description  is  particularly  useful  for  analysis  and  synthesis  of 
hybrid  modulation/demodulation  systems  which  are  designed  to  have  well-defined 
range  responses  and  which  may  be  period-modulated. 
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Figure  18:  System  2 Block  Diagram 
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Figure  19:  System  2 Graphical  Mapping  Diagram 
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Figure  23:  System  3 Measured  Range  Respons 


PH  INPUT 


Vi+  -r  Cj  + 7 C;>+ 


Figure  26  : System  4 Block  Diagram 
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Figure  27  : System  4 Graphical  Mapping  Diagram 
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Figure  30:  CW  Range  Response  Showing  Transmission 
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